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Abstract  Previously,  the  stability  of  a  premixed  flame  in  a  stagnation  flow  was  discussed  for  a  res¬ 
tricted  class  of  disturbances  that  are  self-similar  to  the  basic  undisturbed  flow:  thus,  flame  fronts  with  cchtu- 
gations  only  in  the  cross  stream  direction  were  considered.  Here,  we  consider  a  more  general  class  of  three- 
dimensional  flame  front  perturbations  which  also  permits  corrugations  in  the  streamwise  direction.  It  is 
shown  that,  because  of  the  stretch  experienced  by  the  flame,  the  hydrodynamic  instability  is  limited  only  to 
disturbances  of  short  wavelength.  If  in  addition  diffusion  effects  have  a  stabilizing  influence,  as  would  be 
the  case  for  mixtures  with  Lewis  number  greater  than  one,  a  stretched  flame  could  be  absolutely  stable. 
Instabilities  occur  when  the  Lewis  number  is  below  some  critical  value  less  than  one.  Neutral  stability  boun¬ 
daries  are  presented  in  terms  of  the  Lewis  number,  the  strain  rate  and  the  appropriate  wavenumbers.  Beyond 
the  stability  threshold  the  two-dimensional  self-similar  modes  always  grow  first.  However,  if  disturbances  of 
tong  wavelength  are  excluded,  it  is  possible  for  the  three-dimensional  modes  to  be  the  least  stable  one. 
Accordingly,  the  pattern  that  wilt  be  observed  on  the  flame  fronL  at  the  onset  of  instability,  will  consist  of 
either  ridges  in  the  direction  of  stretch  or  the  more  common  three-dimensional  cellular  structure. 
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ISEBfiOSCZlfiB 

The  response  of  a  premixed  flame  to  stretch  Is  a  problem  of  fundamental 
Importance  in  combustion  theory.  Results  of  studies  on  the  structure, 
propagation  and  extinction  of  stretched  flames  have  been  used  In  experimental 
and  theoretical  Investigations  (Wu  and  Law  1984;  Law  et  al.  1986;  Tien  and 
Matalon  1990;  Dlxon-Lewls  1991)  aimed  at  vinderstandlng  the  structure  and 
burning  characteristics  of  laminar  flames.  Results  of  these  studies  have  been 
also  applied  to  turbulent  combustion,  In  particular.  In  the  "reaction  sheet" 
regime  (Williams,  1985;  Peters,  1986)  where  the  turbulent  flame  Is  regarded  as 
an  ensemble  of  stretched  laminar  flamelets.  Thus,  considerable  attention  has 
been  given  to  stretched  flames  In  the  literature  where  a  large  number  of 
theoretical  and  experimental  papers  can  be  found.  In  spite  of  these  numerous 
Investigations  the  effect  of  stretch  on  flame  front  Instabilities  has  not  been 
completely  clarified. 

Experimental  studies  (Ishlzuka  and  Law,  1982;  Ishlzuka,  1988)  on  the 
stability  limits  of  premixed  flames  clearly  Indicate  the  Importance  of  flame 
stretch.  Observations  were  made  on  propane/alr  premixed  flames  In  an 
axlsymmetrlc  stagnation  point  flow.  The  flow  field  In  such  a  configuration 
Imposes  on  the  bulk  of  the  flame  a  well  defined  strain  rate,  or  stretch, 
controlled  by  the  mass  flow  rate  of  the  Incoming  fresh  mixture.  Depending  on 
the  strain  rate,  various  patterns  appear  on  the  flame  front  for  rich 
mixtures.  For  a  weak  strain  rate  a  three  dimensional  cellular  flame  evolves, 
as  may  be  expected;  cellular  flames  are  commonly  observed  In  experiments 
performed  In  large  tubes  (Clavln,  1985),  neunely  In  the  absence  of  stretch. 

Upon  Increasing  the  strain  rate,  the  flame  approaches  the  stagnation  point  and 
a  different  pattern,  characterized  by  ridges  along  the  radial  direction, 
emerges  on  the  surface  of  the  flame  front.  For  sufficiently  large  values  of 
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the  strain  rate,  all  instabilities  are  suppressed  resulting  with  a  smooth  flat 
flame  surface.  In  contrast,  for  lean  mixtures  the  observations  indicate  an 
absolutely  stable  flame  front  for  equivalent  values  of  the  strain  rate.  These 
experiments  clearly  point  out  the  role  of  flame  stretch,  combined  with 
hydrodynamical  and  dlffusional- thermal  effects,  in  suppressing/promoting 
instabilities . 

A  complete  theoretical  study  on  the  effect  of  stretch  on  flame  stability 
has  not  yet  been  given.  Although  there  is  no  difficulty  in  formulating  the 
linearized  disturbance  equations  for  a  flame  in  a  stagnation -point  flow,  the 
resulting  problem  consists  of  a  coupled  system  of  partial  differential 
equations  with  variable  coefficients.  These  coefficients,  which  depend  on  the 
unperturbed  base  flow,  can  be  written  explicitly  only  for  a  hydrodynamic 
model,  where  the  flame  is  treated  as  a  surface  of  discontinuity  (Eteng  et  al., 
1986).  Furthermore,  since  the  basic  state  is  a  non-parallel  flow,  with 
velocities  varying  normal  and  tangential  to  the  flame  front,  the  disturbance 
equations  are  not  separable.  Thus,  the  linear  stability  problem  cannot  be 
resolved  using  normal  modes  in  both  directions  that  are  parallel  to  the  flame 
front.  An  exception  is  a  class  of  disturbances  self-similar  to  the  basic 
state  which  has  been  considered  by  Kim  and  Matalon  (1990).  In  this  case, 
perturbations  have  periodic  wavlness  in  the  cross -stream  direction  only  with 
no  corrugations  in  the  direction  of  flow  divergence.  Self  similar 
disturbances  were  first  considered  by  Gdrtler  (1955)  and  more  recently  by 
Wilson  and  Gladwell  (1978) ,  in  analyzing  the  stability  of  the  incompressible 
viscous  stagnation-point  flow  without  chemical  reaction.  More  general 
disturbances  were  later  admitted  by  Spalart  (1989)  and  firattkus  and  Davis 
(1990)  in  dealing  with  the  fluid  dynamical  problem.  Spalart  solved  the 
initial  value  problem  numerically  and  found  that  the  long  term  behavior  of  the 
solution  approaches  the  spatial  structure  of  the  self  similar  mode.  Brattkus 
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and  Davis  expressed  the  disturbances  in  terms  of  an  infinite  series  of  coupled 
modes  and  concluded  from  their  analysis  that  the  least  stable  mode  is  that 
corresponding  to  a  self-similar  disturbance.  Although  these  studies  provide 
some  Justification  for  the  use  of  self-similar  disturbances,  perturbations 
with  periodic  waviness  in  the  streamwlse  direction  are  not  of  negligible 
significance,  especially  for  the  combustion  problem.  These  modes  are 
undoubtedly  Important  in  the  nonlinear  development  of  flame  disturbances  and 
must  be  included  in  any  description  of  the  three  dimensional  cellular  pattern 
which  has  been  observed  in  experiments. 

In  addition  to  the  analysis  of  Kim  and  Matalon,  the  only  other  reported 
study  on  the  stability  of  a  flame  in  a  stagnation-point  flow  is  that  of 
Sivashinsky  et  al.  (1982).  Their  analysis,  however,  is  restricted  to  a  very 
weak  stretch  rate.  Furthermore,  they  suppressed  the  hydrodjmamic  instability 
by  neglecting  the  effect  of  thermal  expansion,  thus  deriving  results  that 
pertain  to  the  effect  of  stretch  on  the  diffusional- thermal  instability  only. 

In  the  present  work,  we  analyze  the  stability  of  a  premixed  fl£une  in  a 
plane  stagnation-point  and  consider  a  more  general  class  of  three-dimensional 
disturbances  which,  in  particular,  allow  for  flame  front  corrugations  in  the 
streamwlse  direction.  As  is  well  known,  the  hydrodynamic  instability  (Landau, 
1944;  Darrleus,  1945)  does  not  permit  the  existence  of  premixed  flames  that 
are  too  flat.  Although  diffusion  often  introduces  stabilizing  influences,  as 
is  the  case  for  mixtures  with  Lewis  number  greater  then  one,  long  wavelength 
disturbances  are  hydrodynamically  unstable.  Our  present  aim  is  to  determine 
whether  a  sufficiently  large  stretch  can  stabilize  these  long  wavelength 
disturbances.  We  have  also  incorporated  thermo -diffusive  effects  in  our  model 
in  order  to  determine  the  extent  that  a  flame  needs  to  be  stretched  to  remain 
stable.  For  two-dimensional  self-similar  disturbances,  Kim  and  Matalon  (1990) 
found  that  stretch  indeed  stabilizes  the  flame.  In  particular,  they  provided 
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the  modification  resulting  from  stretch,  to  the  neutral  stability  curve  drawn 
previously  for  freely  propagating  flames  (Matalon  and  Matkowsky,  1984;  Jackson 
and  Kapila,  1984),  thus  identifying  the  critical  conditions  at  the  stability 
threshold.  Here,  we  shall  not  only  extend  their  results,  but  also  identify 
the  effect  that  such  disturbances  have  on  the  pattern  that  evolves  on  the 
flame  front  and  in  the  induced  flow  field  when  an  instability  develops.  It  is 
shown  that  beyond  the  stability  threshold,  the  two-dimensional  self-similar 
modes  always  grow  first.  However,  if  disturbances  of  long  wavelength  are 
excluded,  it  is  possible  for  the  three-dimensional  modes  to  be  the  least 
stable  one.  Accordingly,  the  pattern  that  will  be  observed  on  the  flame  front 
at  the  onset  of  the  instability,  will  consist  of  either  ridges  in  the 
direction  of  stretch,  or  the  more  common  three-dimensional  cellular  structure. 
These  results  appear  in  agreement  with  the  experimental  observations 
summarized  above,  though  a  direct  comparison  requires  analyzing  the 
axlsymmetrlc  configuration.  Having  established  a  procedure  for  testing  self- 
similar  base  flows  for  stability  against  general  disturbances,  we  shall 
discuss  the  stability  of  a  premixed  flame  in  an  axlsymmetrlc  stagnation-point 
flow  in  a  future  publication. 
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FQWWlAnOW 

The  mathematical  formulation  Is  identical  to  that  of  Kim  and  Matalon 
(1990).  A  combustible  mixture  impinges  toward  a  flat  surface  held  at  x  >  0 
(see  Fig.  1).  Far  upstream,  the  flow  field  is  the  classical  stagnation-point 
flow  characterized  by  the  strain  rate  e.  The  flame,  which  separates  the 
burned  products  from  the  fresh  unbumed  mixture,  is  considered  thin  and  is 
therefore  represented  by  the  surface  ♦(x.y.z.t)  -  0,  where  ♦  >  0  is  the  burned 
gas  region.  The  flame  structure  resolved  on  the  scale  of  the  diffusion  length 
Ljj  -  where  D.j.  is  the  mixture  thermal  diffusivity  and  the  laminar 

flame  speed,  yields  jump  conditions  for  mass  and  momentvus  across  the  flame,  as 
well  as  an  equation  describing  the  evolution  of  the  flame  front  (Matalon  and 
Matkowsky,  1983).  Since  viscosity  has  only  a  secondary  effect  on  flame 
stability,  it  will  be  neglected  everywhere.  Consequently,  the  flow  field  on 
either  side  of  the  flame  is  described  to  all  orders  by  the  inviscid, 
incompressible  Euler  equations  but  with  distinct  densities  in  the  unbumed  and 
burned  gas.  The  jump  conditions  and  the  flame  speed  equation  are  only  correct 
to  0(5),  where  5  ■  ^jj/L  «  1  is  the  ratio  of  the  diffusion  length  to  a 

hydrodynamic  length  L. 

In  writing  the  governing  equations  the  length  L  is  used  as  a  unit  of 
distance,  the  laminar  flame  speed  as  a  unit  of  speed,  L/S^  as  a  unit  of 

time  and  the  state  of  the  fresh  mixture  p^,  Tq  is  used  to 

nondimens ionalize  the  density  p,  pressure  p  and  temperature  T,  respectively. 
These  equations  are; 
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for  ♦(x,y,z,t)  <  0 


V*v  -  0 


across  ♦(x,y,2,t)  -  0 

[p(v.n)  .  Vjl  -  fi(fn  a)  K  (3) 

t 

[v  X  n]  -  -  «/9|{Vxv]  +  2(0-1)  VxnJ-  (4) 

[p  +  p(v.^)(v.n  -  V^)l  -  -  (0-1) (1+^)  V.n  +  V^(in  o)  k]>  (5) 

where  I-]  denotes  the  jump  in  the  quantity,  i.e.  (O^j^j-ned  '  ^'^unburned' 


at  ♦(x,y,z,t)  -  O' 


v.n  -  -  1  -  S(0^yJi/2)  K 


In  the  above  D/Dt  is  the  convective  derivative  and  the  variable  p  stands  for 
the  small  pressure  deviations  from  the  ambient  pressure  level  (of  the  order  of 
a  squared  Mach  number) .  The  density  is  given  by 


J  1  ♦  <  0 

\l/a  «  >  0 


where  a  is  the  theirmal  expansion  parameter  (o  >  1).  The  velocity  of  the  flame 
surface  and  the  unit  normal  vector  n  are 


^f  - 


_L_  ai 
I7*i  at  ’ 


|V*| 


The  flame  stretch  is  given  by 


K  -  V.n  -  n.  Vx(vxn) 


(9) 
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with  the  right  hand  side  evaluated  at  •(x,y,z,t)  -  O'.  The  remaining 
coefficients  are 


i  -  (E/RT^XLe-D/a^ 

where  E  is  the  activation  energy  of  the  chemical  reaction  and  R  the  gas 
constant.  Note  that  and  7  are  both  positive  quantities  which  depend  solely 
on  the  thermal  expansion  parameter,  and  i  is  a  scaled  Lewis  number  that 
measures  the  deviation  of  Le  from  one.  Finally  the  boundary  conditions  are 
u  -  0  at  X  -  0  (10) 

V  -  (-€x,€y,0)  as  x  -►  -«>  (11) 


where  v  -  (u,v,w)  and  e  -  (L/S^)e  is  the  dimensionless  strain  rate. 

The  basic  state  corresponds  to  a  steady  planar  flame  front  located  at 
X  -  -d.  The  velocity  field  is  given  by  (Eteng  et  al.,  1986) 

V  -  {  F(x),  -yF'(x),  0  ),  (12) 

where  F(x)  -  Fq(x)  +  5Fj^(x)  and  prime  denotes  differentiation  with  respect  to 

X.  Here, 


Fo(x)  - 


'-c(x  +  Sq) 
-ct(^)^/^  sin  qx 


-CO  <  X  <  -d 
-d  <  X  <  0 


F^(x)  - 


-  ‘*1 

lij. 
2  ^ 


sin  qx  -  qx  cos  qx 


-CO  <  X  <  -d 

-d  <  X  <  0 


where  q 


P 


tja-l/a.  The  pressure  field  takes  the  form 


-CO  <  X  <  d 

-d  <  X  <  0. 


(13) 
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The  flame  standoff  distance  d  and  displacement  distance  a  are  expressed  as 
d  -  dQ  +  5dj^  and  a  -  ag  +  5a^,  where 

do  “  ’?/*  •  dj^  -  -  ^  -  ytn/^ 

Sq  -  dg  -  1/e  ,  aj^  -  dj^  +  ^  +  il/l 

with  ij  •  (a  tan* V«7-l)/yo-l.  Note  that  for  a  very  weak  stretch  (e  «  1)  the 
flame  recedes  back  to  the  nearly  uniform  flow,  far  upstream.  For  0(1)  values 
of  e  the  flame  stands  at  a  finite  distance  from  the  plate,  thus  Justifying  the 
inviscid  approximation  assumed  earlier.  Only  when  c  »  1  does  viscous  effects 
become  important  because  the  flame,  in  this  limit,  is  pushed  into  the  thin 
viscous  boundary  layer  adjacent  to  the  plate. 
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•J  M  I'J ;  J  JA>  *  Lt  ii-': 


Ue  Introduce  small  disturbances,  superimposed  to  the  basic  state  (12)- 
(13),  of  the  form 

'  u(x.y) 
v(x,y) 

w(x,y) 

.  P(x,y) 

and  express  the  perturbed  flame  front  as 

X  -  -d  +  ^(y)e^*'^‘^*'  (15) 

where  k  is  the  wavenumber  associated  with  the  spatial  periodicity  in  the  z 
direction  and  u  is  the  growth  rate.  Note,  that  using  (14) -(15),  a  solution 
that  assumes  normal  modes  in  z  and  t,  but  not  in  y,  has  been  sought. 
Substituting  into  the  governing  equations  (l)-(2)  and  then  linearizing  about 


the  basic  state,  one  obtains 

u  +  V  +•  ikw  -  0  (16a) 

X  y 

F(x^  -  yF'(iy  +  (a,  +  F')a  -  -  p^/p  (16b) 

-  y^'K  +  (w  -  F')<^  -  yF"G  -  -  6  //>  (16c) 

Fw  -  yF'w  +  ww  -  -  ikp/p  (16d) 

X  y 


where,  here  and  thereafter,  subscripts  denote  partial  differentiation. 

In  the  absence  of  vorticity  far  upstream,  and  as  a  result  of  the  density 
being  constant,  the  flow  field  up  to  the  flame  can  be  determined  from  a 

2 

velocity  potential  s!>(x,y,x,t)  which  satisfies  7  ip  -  0.  Writing 

A ,  V  ikz+«t 
<p  -  <p(x,y)  e 

Laplace's  equation  yields 

A  ,  A  ,  2  A 

®  +®  -ka  —  0 

’^xx  ^yy  ^ 


(17) 
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from  which  the  velocity  field  can  be  determined  as 


U  -  ,  V  -  (jp  ,  w  -  ik^ 

*  y 

The  pressure  is  obtained  from  any  one  of  the  components  of  the  momentum 
equation,  giving 


(18) 


P  -  yF'^y  -  (19) 

Thus,  disturbances  in  the  unburned  gas  are  determined  in  terms  of  a  single 
variable  $>.  In  the  burned  gas,  however,  the  flow, is  no  longer  potential 
because  the  vorticity  produced  at  the  flame  is  transported  downstream. 
Consequently  the  stability  problem  reduces  to  solving  (17)  in  the  unburned 
region  and  the  system  (16)  in  the  burned  region. 

Expressing  the  jump  conditions  (3) -(5)  and  the  flame  speed  relation  (6) 
at  the  mean  position  x  -  -d^,  then  linearizing  about  the  basic  state  and  using 


the  solution  expressed  in  terms  of  Ip  for  the  unburned  region,  one  obtains  the 
following  conditions: 

At  X  -  -  dp: 

n  -  -  6(<7-l)<J[^+(l+f7)J^j  +  jfe^»7(y^)y  +  2^[^yy-  (20a) 

V  -  Vy  -  2^n|^(l+a/n)y^  + 

+  «3^yy  +  (‘-H-O^y  +  ^yVy}  (20b) 
w  -  ikv  -  ik(a-l)0  -  iki^^  ^»?'^(w- «<?/»> )^  +  +  uxp  +  (20c) 

P  ”  -ey^y  -  [l-«e(l-»7)^jvjj  +  «(o-l)|j^l+/9+^(»?+2)j  l^^yy  -  k^^fj 

+  T-*[vyy  -  +  ^qj^(2€^-«^)^  -  2ew  y^  -  *^y^^yy]}  (20d) 

+  ey^y  +  *{' (^-"2^^  [^yy  '  +  2^(’^'l>[^yy  '  (2°®) 
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Far  upstream,  l.e.  as  x  -<»,  all  velocities  must  vanish  and  therefore  ^  must 
be  bounded.  At  the  plate,  only  the  axial  velocity  must  vanish,  i.e.  u  -  0  at 
X  -  0,  as  appropriate  for  inviscid  flows. 

The  system  (16)  for  x  >  -d  together  with  equation  (17)  for  x  <  -d,  the 
relations  (20a-e)  at  the  flame  front  and  the  boundary  conditions  at  x  -  -<»  and 
at  X  -  0,  constitute  an  eigenvalue  problem  for  the  determination  of  the  growth 
rate  u.  The  eigenvalues  u  and  the  associated  eigenfunctions  depend  on  the 
following  parameters:  the  thermal  expansion  a,  the  strain  rate  c,  the 
(reduced)  Lewis  number  f,  the  flame  thickness  6  and  the  form  of  the 
disturbance  represented  by  k,  the  wavenumber  associated  with  the  z-direction, 
and  an  appropriate  parameter  that  remains  to  be  determined,  identifying  the 
form  of  the  corrugations  in  the  y- direction. 


As  already  pointed  out  In  the  introduction,  the  eigenvalue  problem  we 


have  just  obtained  is  not  separable  so  that  normal  modes  in  y  are  not 
admissible  solutions.  An  exception  is  the  self  similar  form 


{  u,  V.  w,  p  )  -  {  U(x),  yV(x),  W(x),  P(x)  ) 


(21) 


with  ^  -  const.,  considered  by  Kim  and  Matalon  (1990).  In  this  case,  all 
variables  are  functions  of  the  axial  coordinate  x,  with  the  exception  of  the 
streamwise  velocity  component  which  is  linearly  proportional  to  y.  Clearly, 
this  follows  from  the  fact  that  the  stagnation-point  flow  is  a  local  solution 
valid  near  the  axis  (y  -  0) .  It  is  plausible,  therefore,  to  expect  that  more 
general  disturbances  will  include  higher  powers  of  y  as  well. 

Following  Brattkus  and  Davis  (1990),  further  insight  in  choosing  the 
allowable  disturbances  can  be  obtained  by  examining  the  vorticity  equation 

A  A  A 

uO  +  F  0  -  yF’  n  -  -ikyF"  w  +  F"  v  +  yF'  "  u  ,  (22) 

X  y 


A 


obtained  by  cross-differentiating  equations  (16b)  and  (16c),  where  0  -  v^ 

is  the  z-component  of  the  vorticity  perturbation.  Note  that  the  base  flow 

contains  vorticity  in  the  burned  gas  O  -  -yF",  which  gives  rise  to  the 

s 

nonzero  terms  on  the  right  hand  side  of  equation  (22).  If  one  Ignores  0 

s 


and 


considers  a  potential  flow  F  •  -x,  the  right  hand  side  of  equation  (22) 
vanishes.  The  simplified  equation,  for  -d  <  x  ^  0,  -*  <  y  <  «,  admits 

A  n+w  n 

separable  solutions  of  the  form  0  ~  x  y  .  Depending  on  n,  and  on  whether 
the  disturbances  grow  or  decay,  the  vorticity  may  either  vanish  or  amplify  as 
X  -»  0;  in  the  latter  case  the  vorticity  decays  in  the  boundary  layer  near  the 
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wall  due  to  viscous  action.  In  any  case  O  Increases  with  distance  from  the 
axis  as  y". 

In  the  presence  of  vortlclty  In  the  base  flow,  equation  (21)  Is  no  longer 
separable  and  therefore  the  disturbances  are  not  expressible  In  a  simple  form 
as  In  (22).  Nevertheless,  we  may  write 


A 


U 


00 


n-0 


00  00 
V  -  S  v„(x)y",  w  -  E  lkw„(x)y" 
n-0  "  n-0  " 


fp 


E 

n-0 


»-  J  A„y 
n-0 


n 


(23) 


expecting  that  now  the  various  modes,  corresponding  to  different  Indices  n, 
are  all  coupled.  Brattkus  and  Davis  (1990)  have  also  expressed  the 
disturbances  as  an  Infinite  series,  but  In  terms  of  Hermlte  polynomials  which 
arise  naturally  when  viscous  effects  are  Incorporated  In  the  vortlclty 
equation.  Finally,  It  should  be  pointed  out  that  unlike  the  self* similar 
disturbances  (21) ,  which  are  essentially  two  dimensional  with  corrugations 
only  In  the  z-dlrectlon,  the  disturbances  considered  here  are  three 


dimensional . 


Upon  substituting  (23)  Into  the  governing  equations  (16) -(17)  and  the 


boundary  conditions  (20),  the  following  system  results: 
for  X  <  -d^: 

*n  ■  "  -(n+l)(n+2)*^^2  <24) 

for  -d^  <  X  <  0: 

+  (n+l)V^^^  -  -  0  (25a) 

FU^  +  («  -  (n-l)F*]U^  -  -  oP^  (25b) 

FV^  +  [w  -  (n+l)F']Vjj  -  -  a(n+l)P^^^  +  (25c) 

FW'  +  tw  -  nF']W„  -  -  oP  (25d) 

n  n  n 


Un  -  -  6(o-l){[^+(U,)f][(n+l)(n^2)A^^2  *  ^\]  ^ 

+  J^[(n+l)(n+2)*^^2  *  ^\]}  <26a) 

''n  -  ("^DVl  •  *  '\-l]  ■ 

+  [<^€(n+l)-«a/f;)](n+l)A^^j^  +  [t.H-(n+l)c  ]  (n+l)*^^j^j^  (26b) 

^n  "  *n*  2^«7|[<^<n-ea/i7)JA^  +  ((4H-ne)#^|-  (26c) 

P^  -  -((^n«)*^  -  [^l-6e(l-i7)j2]*;  +  6(a-l)-[7-«[(n+l)(n+2)*^^2  * 

+  [l+^4.|^(i,+2)][(n+l)(n+2)A^^2  * 

+  ^»f[^2€^-w^-2n€w-€^n(n-l)jA^|-  (26d) 
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-  [(.H-e(n+l)]A^  +  fi|-[i8+J^][(n+l)(n+2)A^^2  ' 

+  f^(»J-l)[(tv+l)(n+2)*^^2  *  ^26e) 

at  »  ■  Q 

-  0  (27a) 

at  X  -  -»  -«e 

are  bounded  (27b) 

We  note  that  the  problem  for  the  n^^  order  mode  is  coupled  to  the  systems 
corresponding  to  the  n-1,  n-fl  and  n+2  modes.  The  mode  coupling  in  the 
equations  appears  only  on  the  right  hand  side  of  equation  (25c)  and  is  a 
result  of  downstream  pressure  gradients,  and  of  vorticity  produced  by  the  base 
flow  at  the  flame  front.  Mode  coupling  in  the  conditions  (26)  is  of  two 
origins:  The  first,  which  is  a  result  of  vorticity  generation  at  the 
perturbed  flame  front,  appears  in  the  0(1)  terms  on  the  right  hand  side  of 
equation  (26b).  The  second,  which  is  a  result  of  flame  front  curvature  and 
flame  stretch  (effects  that  are  of  the  same  order  as  the  flame  thickness), 
appears  in  the  0(5)  terms  of  these  relations. 

A  direct  observation  suggests  that  the  eigenfunctions  of  the  system  (24)- 
(27)  fall  into  two  distinct  classes,  with  the  velocity  v  being  either 
symmetric  or  anti -symmetric  about  the  axis  y  -  0.  Each  of  these  two  classes 
has  its  own  set  of  eigenvalues.  We  refer  to  the  even  or  odd  solution  as  the 

series  containing  the  sequence  |  ^2^’  '^2n’  ^2n-l’  '^2n  |  "“0*1*2...  ^or* 

{  ^2n-l*  ^2n-l’  '^2n'  '^2n-l  |  "  “  }  ,  respectively.  The  classification 

is  made  according  to  whether  the  disturbed  flame  front  ^  and  hence  u  and  ^  are 
even  or  odd  functions  of  y;  v  always  has  the  opposite  parity.  We  determine 
these  solutions  by  truncating  the  series  at  n  -  N  and  finding  the  eigenvalues 
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at  any  given  level  of  truncation.  For  a  given  N,  one  of  the  two  following 
system  of  equations  results  for  the  determination  of  the  eigenvalues: 

1 


^N-1^ 

2 

NV  -  k  U  — 

”  N  N-1 

0 

(28a) 

FU'.i  +  [«  - 

(N-2)F']Ujj  ^  - 

-  "*’n-i 

(28b) 

FV^  +  [« 

-  (N+l)F']Vj^  - 

^"«N-1 

• 

(28c) 

FW' .1  +  [CO  - 

(N-l)F'lWj,.^  - 

■  "'h-i 

(28d) 

System  II 


FU^  +  [«  -  (N-l)F']Uj,  -  - 
FV'.i  +  [«  .  NF'JV^.^  -  -  +  F"Uj,.2 

FW^  +  [w  -  NF' ]V^  -  -  aP^ 


(29a) 

(29b) 

(29c) 

(29d) 


The  remaining  equations,  for  the  variables  with  lower  indices,  are  only  needed 
for  the  determination  of  the  eigenfunctions.  We  note  parenthetically  that 
equation  (29c)  is  uncoupled  from  the  remaining  equations  in  system  II  and  is 
therefore  not  needed  for  the  determination  of  w.  Each  of  these  two  systems  is 
to  be  supplemented  by  the  conditions  (26a-d),  truncated  appropriately.  The 
latter  requires  the  determination  of  4^^  and  however,  from  equation  (24) 

one  notes  that  the  solution  that  decays  as  x  -»  •<»  always  satisfies  4^  -  k4^ 

which,  when  substituted  in  (26e)  determines  4„. 

As  N  increases  from  N  -  0,  we  first  find  the  solution  Vq  -  0;  for  N  -  1 
we  find  the  first  even  mode  (Uq;  V^;  W^)  which  is  the  self  similar  solution 


(see  eqiiatlon  21)  and  the  first  odd  mode  (Uj^;  VqI  Wj^)  ;  for  N  -  2  we  find  the 

second  odd  and  even  modes  (Uj^;  V2,  Vq;  and  (Uj,  Uq;  Vj^;  W2,  Wq) 

respectively,  etc. . .  In  general,  when  N  is  odd  system  I  determines  the 
eigenvalues  of  the  N-even  mode  and  system  II  determines  the  eigenvalues  for 
the  N-odd  mode;  when  N  is  even  system  I  determines  the  eigenvalues  of  the  N- 
odd  mode  and  system  II  determines  the  eigenvalues  for  the  N-even  mode. 
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STABILITY  OF  A  STKETCHED  FLMg 

The  mathematical  formulation  presented  above  is  valid  for  j  «  1,  so  that 
the  solution  must  be  expressed  in  power  series  of  S .  Here  we  only  follow  the 
first  two  terms  in  this  expansion  and,  in  particular,  express  the  dispersion 
relation  as 

«  -  UQ(o,e;k,H)  +  « ,i,*k,N)  (30) 

For  a  given  N,  the  eigenvalues  for  the  even  and  odd  modes  are  obtained  by 
solving  equations  (28)  or  (29)  with  their  appropriate  boundary  conditions. 

The  leading  problem,  obtained  by  setting  £  -  0,  yields  The  next  term 

is  obtained  as  a  solvability  condition  from  a  system  with  the  same 
differential  operator  as  that  of  the  leading  problem  but  with  nonhomogeneous 
terms  that  depend  on  the  leading  order  solution.  The  solutions  presented 
below  were  obtained  by  numerically  solving  these  systems  using  the  boundary 
value  solver  COLSYS  (Ascher  et  al.,  1981). 

The  dependence  of  the  solution  on  the  strain  rate  c  can  be  scaled  out  in 
the  whole  problem,  as  shown  by  Kim  and  Matalon  (1990).  Consequently,  the 
dispersion  relation  (30)  can  be  rewritten  as 

w  -  eaQ(<7;K,N)  +  fi«^aj^(a,i;ic,N)  (31) 

where  le  -  k/< . 

The  leading  order  term  in  the  dispersion  relation  corresponds  to  the 
modification  of  the  Darrieus- Landau  instability  resulting  from  flame  stretch. 
For  very  small  e,  the  flame  is  located  far  upstream  in  a  nearly  uniform  flow. 
Thus,  the  hydrodynamic  disturbances  grow  in  the  infinitely  long  downstream 
region  of  burned  gas  according  to 


(32) 
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Landau,  1944;  Darrleus,  1945).  For  moderate  values  of  c,  the  flame  is  located 
at  the  mean  position  x  -  -d.  Consequently  the  hydrodynamic  disturbances  which 
are  now  confined  to  the  region  -d  <  x  <  0  of  burned  gas,  must  vanish  at  a 
finite  distance  from  the  flame  front.  This  damping  effect  is  one  source  of 
stabilization,  the  only  one  experienced  by  the  two  dimensional  self- similar 
disturbances.  An  additional  source  of  stabilization  comes  from  the  stretching 
that  three  dimensional  disturbances  experience  in  the  direction  of  flow 
divergence.  These  two  effects  modify  the  growth  rate  (32),  as  shown  in 
Fig.  2,  where  the  leading  order  term  of  the  dispersion  relation  is  plotted 
against  k  for  a  selected  number  of  even  and  odd  modes  corresponding  to  various 
N.  [Here,  and  in  all  following  figures,  solid  curves  result  from  solutions  of 
system  1  while  dashed  curves  result  from  solutions  of  system  II;  the  notation 
associated  with  the  different  modes  consists  of  a  digit  which  stands  for  the 
value  N  and  the  letters  E  for  even  and  0  for  odd] .  We  note  from  Fig.  2a  that 
disturbances  with  long  wavelength  in  the  z-directlon  (i.e.  with  small  k)  are 
now  damped  out.  However,  for  a  given  k,  disturbances  with  short  wavelength  in 
the  y-direction  (i.e.  corresponding  to  large  N)  are  the  more  stable  one.  For 

a  fixed  c  there  exists  therefore  a  value  k~(N)  such  that  the  only  growing 

modes  are  those  with  k  >  k  (N) ,  with  k'  increasing  with  N.  For  comparison,  we 
have  plotted  in  Fig.  2b  the  relative  growth  rate  u>q  compared  to  that  of  an 

unstretched  flame  .  For  a  fixed  k  this  graph  identifies  the  dependence  of 

the  growth  rate  on  1/e .  Note  that  the  ratio  is  less  then  one  for 

e  >  0  and  approaches  one  as  e  -»  0. 

The  effects  of  diffusion  appear  in  the  0(j)  term  of  the  dispersion 
relation  (30).  For  an  unstretched  flame,  the  correction  to  the  growth  rate 


is  of  the  form 
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<7(c-l)(l+^)  +  2a/9(l+ajs)-i  rl  a(l+a«)  (a+a^x)-]!  « 

) - H  ■]}  <“> 

(Pelc6  and  Clavln,  1982;  Hatalon  and  Matkowsky,  1983).  In  the  presence  of 
stretch  It  is  clear  from  the  boundary  conditions  (26)  that  the  solution  also 
depends  linearly  on  i  so  that  one  can  write 

t»l  “  aj^j^(a:ic,N)  +  (32) 

Our  calculations  indicate  that,  similar  to  the  case  c  -  0,  the  coefficients 
and  l^oth  negative  (see  Fig.  3).  This  implies  that  diffusion  has 

always  a  stabilizing  influence  when  the  Lewis  ntimber  Le  is  greater  than  one 
(i.e.  i  >  0) .  Diffusion  is  a  source  of  instability  only  when  Le  <  1;  the 
precise  determination  of  the  critical  Jt  will  be  discussed  later.  Note  that 
the  coefficients  and  quite  insensitive  to  variations  in  N;  for  N 

as  large  as  5  the  variations  in  remain  indistinguishable  while  those  in 

a^2  less  than  10%. 

The  combined  effects  of  hydrodynamic,  including  stretch,  and  diffusion 
can  be  understood  by  considering  the  expansion  (30)  as  an  approximation  to  the 
exact  dispersion  relation.  This  procedure  was  found  to  provide  qualitatively 
correct  results  for  the  case  «  -  0  (Matalon  and  Matkowsky,  1984)  which  were 
also  quantitatively  in  agreement  with  the  numerical  results  of  Jackson  and 
Kaplla  (1984)  for  small  k.  The  results  for  a  stretched  flame  are  shown  in 
Fig.  4  where  the  growth  rate  u>  is  plotted  versus  k  for  a  selected  number  of 
modes  N.  The  stabilization  effect  of  diffusion  is  clearly  seen  when  compared 
to  Fig.  2.  First,  disturbances  with  short  wavelength  in  the  z-dlrectlon  are 

damped  out  so  that,  for  a  given  «,  growing  modes  are  limited  to  k*  <  k  <  k^. 

Second,  both  k  and  k^  decrease  with  increasing  N  so  that,  for  sufficiently 
large  N,  all  modes  are  damped  out  (i.e.  u>  <  0  for  all  k  and  N)  resulting  with 
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a  stable  flame  front.  We  note  from  Fig.  4a,  which  corresponds  to  i  -  -  0.53, 
that  the  most  dangerous  disturbance  is  IE,  namely  the  two-dimensional  self¬ 
similar  disturbance,  with  0.5  <  k  <  3.  However,  for  i  -  -  0.54  (see  Fig.  4b), 
three  dimensional  disturbances  with  N  <  2  and  relatively  much  larger  k,  grow 
at  a  faster  rate  than  mode  IE  does.  As  i  decreases  further,  more  of  the 
three-dimensional  disturbances  grow  faster  than  IE;  for  example,  when 
i  -  -  0.6  all  five  modes  shown  in  Fig.  4b  have  growth  rates  that  are  larger 
than  that  of  IE.  The  implication  of  this  observation  will  be  discussed  next. 

The  neutral  stability  boundary  is  obtained  by  setting  w  -  0  in  the 
dispersion  relation  (31).  The  results,  plotted  in  Fig.  5,  are  presented  in 
the  K-i  plane  for  several  values  of  the  parameter  eS .  The  various  curves 
correspond  to  different  values  of  N.  The  region  to  the  right  of  these  curves 
is  where  u  is  negative  and  the  flame  is  stable.  The  stability  boundary  is  the 
envelope  of  these  curves  and  the  critical  is  defined  as  the  maximum  value 

of  t  taken  over  all  k  and  N.  Note  that  i  <  0  so  that  a  stretched  flame  in  a 

c 

premixed  combustible  mixture  with  Lewis  number  Le  >  1,  such  as  a  lean 
hydrocarbon/air  mixture,  is  stable.  Upon  decreasing  i,  for  example  by  making 
the  mixture  richer,  an  instability  develops.  From  Fig.  5  one  sees  that  at  the 
stability  threshold  the  two-dimensional  self-similar  disturbances  grow  first 
and  that  the  transition  occurs  at  a  larger  value  of  I,  the  smaller  eS  is.  The 
pattern  that  will  be  observed  on  the  flame  front,  in  this  case,  is  one  with 
ridges  in  the  z- direction  only.  If  the  finite  width  of  the  apparatus  in  the 
laboratory  exclude  disturbances  with  long  wavelength,  it  is  possible  for  a 
three  dimensional  mode  to  be  the  least  stable  one  and,  consequently,  for  a 
cellular  flame  to  be  observed  at  the  onset  of  instability.  For  example,  if 
disturbances  with  ic  <  3  and  N  <  N°  are  excluded  for  ej  -  2,  one  of  the  dashed 

curves  characterized  by  some  k  and  N  >  N**  will  be  crossed  first  upon 
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decreasing  i.  Thus,  depending  on  the  strain  rate,  two  or  three  dimensional 
patterns  could  result  on  the  flame  front  of  rich  hydrocarbon/air  mixtures. 

The  structure  of  the  flow  field  at  the  onset  of  instability  is  quite  different 
in  these  two  cases,  as  shown  in  Fig.  6.  In  the  figure,  streamlines  of  the 
perturbed  flow  field  on  either  side  of  the  flame  are  plotted  in  the  y  plane. 
Case  (a)  corresponds  to  the  IE  solution  derived  from  system  I  (eqs.  28) 
associated  with  the  right  most  point  of  the  solid  curve  in  Fig.  5;  case  (b) 
corresponds  to  the  10  solution  derived  from  system  II  (eqs.  29)  associated 
with  the  right  most  point  of  the  dashed  curve  in  Fig.  5.  We  note  that  in  the 
latter  case  the  flow  pattern  consists  of  several  arrays  of  vortices  stacked  in 
the  upstream  direction  and  periodically  distributed  in  the  spanwise  direction 
2.  This  whole  structure  is  convected  along  the  diverging  streamlines  of  the 
mean  flow. 
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CQHCLPSIOilS 

We  have  considered  the  stability  of  a  premixed  flame  In  a  stagnation- 
point  flow  against  a  class  of  three  dimensional  flame  front  perturbations.  It 
has  been  concluded  that  the  least  stable  mode  corresponds  to  disturbances 
which  are  self- similar  to  the  basic  state,  in  accord  with  the  known  stability 
results  for  a  viscous  incompressible  stagnation -point  flow  without  chemical 
reactions.  If  disturbances  of  long  wavelength  are  excluded,  however,  it  is 
possible  for  a  three- dimens loanl  mode  to  grow  first.  This  of  course 
determines  the  pattern  that  is  likely  to  be  observed  on  the  flame  front  at  the 
onset  of  instability.  We  have  also  determine  in  this  study  the  stability 
criteria  in  terms  of  the  strain  rate,  the  Lewis  number  and  the  thermal 
expansion  parameter. 
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Fig.  1:  Schematic  representation  of  a  flame  In  a  stagnation-point  flow. 


The  modification  of  the  Darrieus - Landau  instability  resulting  from 
flame  stretch:  (a)  the  growth  rate  for  a  stretched  flame;  (b)  the 
relative  growth  rate  compared  to  an  unstretched  flame.  The  solid 
curves  correspond  to  solutions  of  system  I  with  the  mode 
classification  given  in  the  second  column;  the  dashed  curves 
correspond  to  solutions  of  system  II  with  the  mode  classification 
given  in  the  first  column; 


The  coefficients  and  0^2  calculated  for  even  and  odd  nodes  with 
N  -  1,2, 3, 4, 5  and  <7-5.  Note  that  is  practically  independent 
of  N  and  that  0^2  varies  sli^tly  with  N  depending  on  whether  the 
eigenvalues  result  fron  systen  1  or  II. 


Streanlines  of  the  perturbed  flow  field,  on  either  side  of  the  flane 
located  at  x/d  -  *1,  at  the  onset  of  instability:  (a)  from  the 
solution  IE  derived  from  system  I  for  o  -  5,  <6  >0.8,  i>  -  0.43 
and  K  -  0.8;  (b)  from  the  solution  10  derived  from  system  II  for 

0-5,  e«  -  0.8,  i  -  -  0.53  and  a  -  4.4. 
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